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Abstract We study the relationship between mixed stochastic differential equations and the corre¬ 
sponding rough path equations driven by standard Brownian motion and fractional Brownian motion 
with Hurst parameter H > 1/2. We establish a correction formula, which relates both types of equa¬ 
tions, analogously to the Ito-Stratonovich correction formula. This correction formula allows to 
transfer properties, which are established for one type of equation to the other, and we will illustrate 
this by considering numerical methods for mixed and rough SDEs. 
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1 Introduction 

In this manuscript we will consider two types of stochastic differential equations (SDEs), the so- 
called mixed SDEs (see e.g. mm) and rough SDEs (see e.g. mm driven by standard Brownian 
motion and fractional Brownian motion with Hurst parameter H > 1/2. The mixed SDE reads as 

X^^ = xo+ [ a^{Xf)ds+ [ b^{Xf)d'Ws+ [ c^{Xf)dB^, t>0, (1) 

Jo Jo Jo 

while the corresponding rough path equation is given by 

X^ = xo+ [\^{X^)ds+ [\^{X^)d^Ws+ fc^{X^)d^B^, t>0, (2) 

Jo Jo Jo 

where W = {Wt)t>o is an m-dimensional Wiener process, is an Tdimensional 

fractional Brownian motion with Hurst index H > 1/2, xq , and the coefficients , a^: —>■ 
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b^,b^: —>■ ,c^: R'^ R*^^^ satisfy suitable smoothness assumptions. The precise 

definitions are given below. 

The difference between both equations is the definition of the stochastic integrals with respect to 
the Brownian motion. When dealing with mixed equations, Jq b^{X^)d^Ws is understood as ltd 
integral, while for rough equations b^(X^)d^Ws corresponds to a Stratonovich integral. Both 
equations have been well studied so far; in the unique solvability of the mixed SDE was shown, 
while the well definedness and unique solvability of the rough path SDE have been obtained e.g. in 


m 

In this manuscript we establish a correction formula (see Sections [3] and 0]) between equations 
dll) and ([2]), which extends the Ito-Stratonovich correction formula that goes back to two articles 
of W. Wong and M. Zakai ([U H]). This allows to transfer results valid for rough SDEs to mixed 
SDEs and vice versa. We will illustrate this by establishing the smoothness of the solution map and 
by recovering a limit theorem for the mixed SDE ([T|) in Section [5l furthermore, we point out how 
this limit theorem can be used to construct and analyse numerical methods for mixed SDEs, and 
we show that the “natural” Euler scheme for the rough SDE converges to the rough solution. 


2 Preliminaries 


2.1 Notation and Definitions 


In what follows we will work on a filtered probability space (D, X, (Jd)t>0) IP)) which is rich enough to 
contain all the objects defined below. Let W = {W^^\ ..., be a standard m-dimensional 

Wiener process, ..., be an ^-dimensional fractional Brownian motion 

(fBm) with Hurst index H E (1/2,1), that is a collection of centered, independent Gaussian pro¬ 
cesses, independent of W as well, with covariance function 

Rnit, s) = ^ -\t- , s, t > 0. 

The Kolmogorov theorem entails that fBm has a modification with y-Holder sample paths for any 
'y < H and we will identify B^ with this modification in the following. 

We will use the following standard notation: | • | stands for an absolute value of a real number, the 
Euclidean norm of a finite dimensional vector or of a matrix. Moreover, for a function /: [a, 6] —>■ R 
we define the following (semi-)norms: 


2,b] = sup |/(x)|, 

x£[a,b] 


7, [a,6] - sup 

yG[a 
xjty 


\f{x)-f{y)\ 

x,yG[a,b] 1^ 


7,00,[a,6] 


7,[a,6] 


+ 


If there is no ambiguity we will omit an interval index [a,b]. Eor a vector valued function / = 
(/i,..., /rf): [a, 6] —>■ R'^ a corresponding (semi-) norm is defined as the sum of the (semi-) norms of 
the coordinates /j. Next, for a function /: [a, 6]^ —>■ R, that vanishes on a diagonal, i.e. = 0 

for t E [a, 5], we set 


l7,[a,b]2 


sup 

i,s€[a,6] 

t^s 


It-sD' 
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The set of such functions with a hnite norm ||/||^ ^2 is denoted by C 2 ([a, that is C 2 ([a, 6]^) = 

{/: [a,6]2 ^ M | = 0, t € [a,b], ||/||^j„_fe ]2 < 00 }. 

Moreover, we will use the notation for functions / : —)■ which are 

bounded, /c-times differentiable with bounded derivatives and whose fc-th derivative is Holder con¬ 
tinuous of order 5 > 0. Finally, if a map c : ^ is fixed, we introduce the differential 

operators = Yli=i<^ifii-)dxr 

2.2 Integrating with respect to standard Brownian motion and fBm with H > 

1/2 

For basic facts for Ito or Stratonovich integration with respect to standard Brownian motion we 
refer e.g. to [TOl ITT] . 

The integral with respect to a fractional Brownian motion with Hurst parameter H > 1/2 is 
understood in the pathwise Young sense. Namely, for / € C'^([a,6];R) and g € C^([a,6];R) with 
i' + g, > 1 the integral f{x)dg{x) can be defined as the limit of its Riemann sums and satishes 
the so-called Young inequality 

rb 


J a 


ifis) - fia))dg{s) 




y,[o,r]llfl'llAt,[o,T]l& - 


for all a,b & where > 0 is a constant independent of / and g. Thus, the integral 

f{s)dB^’^^^ for a function / G C^([a, 5]; R) is well defined provided that /3 > 1 — H. More details 
can be found e.g. in mm]- 


2.3 Mixed SDEs 

The mixed equation ([T]) reads as 


X^ = 


pt ^ pt 

X0+ a^(Xf)ds + V / )d'lTi^) 

Jo , 1 Jo 


(3) 


i=i 

i 




(xf )dHf’(^3 ^ t G [0, T], 


where denotes the j-th. column of a matrix. As mentioned previously the integrals with respect 
to the Brownian motions are understood as Ito integrals, while the integrals with respect to the 
fractional Brownian motions are understood as Young integrals. This equation has been analysed 
in a series of articles (HIISIISIE]). The most general result on the existence of a unique solution can 
be found in [6]: 

Theorem 1. Assume that 

(i) G CHR'^;R‘^), i = 1,... ,m, j = 1,... 

(a) i = 1,... ,m, j = 1,... ,i, satisfy a linear growth condition, i.e. there exists 

C > 0 such that 


+ X] ^ \c^J^\x)\ < C(1 + |x|), 

1=1 


X G 


2=1 
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(Hi) sup^gRd |(Dc^^’(-?))(x)| < oo. 

Then equation © has a unique solution, i.e. there exists a unique continuous and (-^t)tg[o,r] adapted 
process X = {Xt)t£[o^T]! which satisfies equation ([T]) for almost all a; E 0. 

The above solution is in fact obtained as the limit (in probability) of the solutions of Ito SDEs 
with random coefficients, namely of 

X^''^ = xo+[ {a^iX^^^) + c^{X^'^)B^’^)ds+ [ b^{X^)d'Ws, tE[0,r], 

Jo Jo 

where = n f*f_i/n)vo ^ ^ [0; ^]) ^ = 1)2,..., is a smoothed fBm, see [5]. 


2.4 Rough paths 

Here we briefly recall some notions of the rough path theory, following the algebraic integration 
approach given in m and the recent monograph |14) . For a detailed exposition the reader is sent 

to [aiattSKii]. 

Definition 1. Let 7 > 1/3. A pair (x,x) = (x^, Xs^t)o<t,s<r € C'>'([0, T]; M”^) x € 2 "^ 
is called a '-y-rough path if = 0 for all t E [0, T] and for aU0<s<u<t<T we have 

Xs,t - Xs,n - = (Xs - Xu) 0 (xt - Xu). 

The function (xs^t)s,t,G[o,T] is called Levy area. 

Remark 1. (ij A Levy area for the fractional Brownian motion B^ with Hurst parameter H > 
1/2 is defined as the collection of Young integrals 

B = (B,,i)o<.<t<T = 1 < i,j < i 


1 0<s<t<T 


(a) The Levy area for standard Brownian can he constructed using Stratonovich integration, i. e. 


W = (yV s,t)o<s<t<T = 


t pu 



odlF® 1 < i,j < mj’ 


0<s<t<T 


where f odW denotes the Stratonovich integral. 


(in) In the same way a Levy area for (t, H/, Ri)t6[o,r] aan be constructed, i.e. using W for the 
iterated integrals of with respect to i,j = l,...,m, and Young integrals for all 

other iterated integrals. 


Together with the notion of the Levy area the following concept is at the core of the algebraic 
integration approach of M. Gubinelli. 

Definition 2. ITe say that a path y E C^([0,r];M^) with v E ( 1 / 8 , 7 ] *•5 a weakly controlled path 
based on x & C"’'([0, Tj; M"*) if the following decomposition holds 


yt-ys = Zs{xt - Xs) + rs,t, 0 < s <t <T, 
with z E C’^{[0,T];R^^^) and r e C|^([0, r]^; ]R*=). 


(4) 
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When (y, z) is a weakly controlled path, then the rough integral of y along x can be defined as 



lim 




— X 




i=i 


(5) 


for i = 1,..., A:, j = 1,..., m, see e.g. Corollary 2 in |13| . Here the limit is taken over all partitions 
P = {0 = t-i = to < ti <...< tn = tn+i = t} such that \V\ = supj^gp \tk — tk-i\ 0. Weakly 
controlled paths are stable under smooth transformations: 

Proposition 1. Let {y,z) be a weakly controlled path based on x with decomposition 0, and let 
(p G C'|(M^;]R”). Then p{y) is a weakly controlled path based on x with decomposition 


Tivt) - pivs) = Zsixt - Xs) + fs,U 


with 

Zs = (D(p){ys)zs, sG[0,T]. 

Using appropriate estimates for the integrals the solution to the rough paths equation 

dyt = (T{yt)dxt, t€[0,r], yo = aeR'^, 

with (T : ^ ]^dxm is obtained via a fixed point argument. 

Theorem 2. Suppose that k € (1/3,7), x: [0;^] —>■ R”^ is a'y-rough path and let a ^ 
such that (2 + 5)'y > 1. Then the equation 

yt = a+[ a{ys)dxs, t G [0,T], 

Jo 


possesses a unique solution in the space of the functions z G C'^([0, T]; R'^) with zq = a. Moreover, 
{y,a{y)) is a weakly controlled path based on x. 

As a consequence of this Theorem and Proposition [1] we have 



= lim 
\p\^o ^ 










i=l 



( 6 ) 


The solution map for a rough equation is locally Lipschitz continuous with respect to the initial 
value and the driving signal. More precisely, we have: 

Theorem 3. Let n G (1/3,7), € C'^’^(R'^;R*^^”^) such that (2 + 5)7 > letx, x be 

'y-rough paths with corresponding Levy areas x, x. Finally, let iyt)t&[o,T]! iyt)te[o,T] solutions 

of the RDEs 


yt = a + 


a{ys)dxs 



tG [o,r]. 
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Then there exist an increasing function Ct'- [ 0 , oo) —>■ [ 0 , oo) such that 


-y|l7,oo,[o,r] < Ct{\\x\\ 


^Il7,oo,[0,r] + Pll7,oo,[0,T] + 11^117,[0,r] + 


X (la — a| + llx — x\ 


ll7,[0,T]) 
7 , 00 , [0,T] + 11^ “ x||2,yj0,T])- 


Returning to our original rough SDE, i.e. to 

rt ^ pt 

X^ = xo+ a^(Xf)ds + V/ (7) 

io Jo 

+ E f t€ [o,r], 

,= 1^0 

the previous results yield a unique solution, if G M*^), i = 1, ... ,m, j = 

for J > 0 arbitrarily small, where the Levy area for {t,Wt, Bt)t£io,T] i® constrncted as in 
Remark [TJ 

3 Prom the rough paths equation to the mixed equation 

Here we show that the solution of the rough SDE ([7]) is the solution of the mixed eqnation ([3]) with 

1 

a^{x) = a^{x) + - b^{x) = b^{x), c^(x) = c^(x), x G (8) 

i=l 

Theorem 4. Let 5 > 0 and c^J^^ G M'^), i = j = Then 

the solution of the rough equation © and the solution X^ of the mixed equation ([3]) with 
coefficients given by ([8]) coincide P-almost surely, i.e. we have 

p(xf = x/^, tG [o,r]) = 1. 


Proof. For the m P £ + 1 dimensional rongh path g = (id, W, B^) denote its Levy area by G = 
{Gs,t)o<s<t<T- Now fix t G [0, T]. Using relation ([6]) for cr = (a^, b^, c ^): ^ x x 

we can write 



lim 




(d 

t^k+l 


l+m+^ 


for all z = 1,..., 1 + m + £. 

Since g^^^ = id and the integrand is continuons, we have 


lim 


E 

ik&'P 


a 


( 1 ) 


i^p)9 


( 1 ) 

^k 5^fc + l 


P—a.s. 



S. 
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For i = 1,... , m, we have 




by definition of the Ito integral since s G [0,r], is bounded and adapted. Moreover, the 

sample paths of are 7 -Holder continuous for all 7 < 1/2 and are Holder continuous of all 
orders X < H, thus we have 

for i = 1,..., by definition of the Young integral. 

Now consider the summands involving the Levy area terms. If {i,j) 0 (2,..., m -|- 1)^, then the 
Young inequality gives 


- 9s'’)<i9u^ 


sup - 

t,5G[0,T] |t - s|l+^ 

t^s 


<00 P — a.s. 


and hence it follows 


lim ^ 0. 

heP 


Next suppose (i,j) € (2,... ,m + 1)^ and i / j, then Gt^^^tk+iihj) = “ Ws^^)d'Wu \ since 

the Stratonovich and the Ito integral coincide due to the independence of IF^*\ Exploiting 

the independence of Gtf.^tk+i from Pt^ obtain 


tkGP 

/‘^/c + 1 

= e[( / (W« - <))dIYi* 

tkGP 

= \Y. E[pi^'V«(Y«)|2]|7:>|2. 

tk^V 

The last term clearly vanishes for \V\ 0. Finally, if f = j we have 

ptk + 1 1 

/ - W,f) o dIF« = -(IT® ^ - IT®)2 


and we obtain 


hm P«u®(Y«)(IT® ^ riP®a®(Yf)ds. 
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The latter follows from 


sup 

sG[0,T] 


n—1 


fc =0 



P—a.s. 


0 


as |P| —0 and a density argument. 

By passing to a subsequence, we deduce that 


m+£+l 


E 

i=l 


(7 


(i) 


pt ^ m 

^ =“-^- ^ {Xf))ds 

m .t I ■ 

+ V / + V 

i=l -^0 i=l 




for all t € [0,r]. Since both sides are continuous in t for almost all a; G the exceptional set can 
be chosen independently of t € [0,7"]. Hence the assertion follows. □ 


4 Prom the mixed equation to the rough paths equation 

Throughout this section we assume that ,b^ ,c^ € C^. Under this assumption the solution 
X^ = to (El) exists, is unique and satisfies E||X||g < oo for all p > 1 and 6 < 1 / 2 , see 

[6]. Now, we will show that X^ is the solution to El with the coefficients 

. m 

a^{x) = a^{x) - b^{x) = b^{x), c^{x) = c^{x), x € (9) 

i=l 


fn order to do to this, we have to show that (X^{(jj),a{X^(uj)) is a weakly controlled path based 
on {g{t){uj),Gt^si^)}o<s<t<T for almost all cj € O, where a = {a^,b^,c^). However, this is a 
consequence of the following two Lemmata. 

Lemma 1. Let h G Then for all 7 G (0,1/2) there exist almost surely finite random 

variables iLr,h, 7 ,g such that 


i: 


{h{X^)-hiXf))dgiu) 


M\ 


< KT,h,i,g ■ |t “ s 


127 


s, t G [0,T], 


where g G {id, with i G (1,... , £}. 

Proof. Since X^ G C'’^([0, Tj) for all 7 < 1/2, the assertion is a direct consequence of the Young 
inequality, which gives 


{h{X:f)-h{Xf))dg{u 




< C„.|* - sr+’'||/{X'“)||,,|„,T,||9||yJ„,T| 


M\ 


for 1 — 7 < 7 ' < 77. 


□ 








Lemma 2. Suppose {Z{t),J^t)tG[o,T] ® stochastic process with 9-Hdlder trajectories for all 6 € 

(0,1/2), such that E||Z||g < oo for all p > 1. Then, for all rj € (0,0), there exists an almost surely 
finite random variable such that 



{Z{v)-Z{s))6' 


W, 


<KT,^lt-sl^/^+^, 


s,t £ [0,T], 


Proof. Let 9 G (0,1/2). Applying the Rodemich-Garcia-Rumsey inequality, we obtain that 


t^S 


f^{Z{v) - Z{s))d'W, 


<Ce, 


v,p 


|a;_y|(l+2r,)p+2 





Now put 

KT,ri = 

The Burkholder-Davis-Gundy inequality gives 


1 . 1 . |i-j,|(i+2i)r+2 ""q ■ 


2p 


{Z{v) - Z{s))d'W, < Ce,p\y - 

Choosing p such that 2(0 — r])p > 1 we obtain 

.r.2p f^E\jy{Z{v)-Z{s))d'W,\^P ^ /■*/■*, . 

- |x-j/|d+2^)p+2 - dxdy<Ce,pJ^J^\x-y\ ^ dxdy < oo. 


□ 


Now we can exploit the representation ([5]) and Proposition [T] to work backwards through the 
proof of Theorem HI which gives: 

Theorem 5. Let G C'|([0,T];R'^), i = l,...,m, j = and suppose that 

= (-T/^)ie[o,T] *-5 the solution to Q, then is a solution to the rough equation ([7]) with 
coefficients given by ([^. 

Note that the smoothness of the drift coefficient of the arising rough path SDE is C/. Within 
the algebraic integration framework it is not known (up to the best of our knowledge) whether such 
an equation has a unique solution. 


5 Application to numerical methods 

5.1 Limit theorem for mixed equations 

At the core of the theory of mixed equations is a limit theorem in [5], which we will briefly recall 
here. Let n G N and define 

rt 

= n B^ds, t >0. 

J {t—l/n)\/0 
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Note that is an (J^i)t>o-adapted Gaussian process such that its trajectories are a.s. differen¬ 

tiable and 


Suppose that = 


(^M,n 


)tg[o,T] ^ solution to 


t > 0. 


vM,n , 

X, ’ = Xo + 


f (a(Xf ■- 

Jo 


) + c(Xf 


)ds+ I 

Jo 


bix^)d'Ws, te[o,T]. 


( 10 ) 


Then we have 


^M,n , re —>■ oo uniformly in probability. (11) 

Setting = (id, W, we can apply the stability of rough paths equations and the relation 

between mixed and rough SDEs to recover and strengthen this result. 

First note that the integral f* c(X^’'^}d^B^’^ coincides with the ordinary Young integral f* c(xf^’’^)B^’’^ds. 
Proceeding analogously to the proof of Theorem |4] we have: 

Proposition 2. Let 6 > 0 and E C'^’‘^(M‘^;M'^), i = 1,... ,m, j = Then the 

solution of the rough equation 

xf’" = xo+ [ a{X^’^)ds+ [ b{X^’^)d^Ws + [ c{X^’^)d^B^’^, tE[0,r], (12) 

Jo Jo Jo 

with 

^ m 

a{x) = a{x) — -'^^'Djj^b^J{x), x E 
i=l 

and the solution of equation (HOD coineide P-almost surely. 

Our aim is now to prove: 

Proposition 3. Let d, b^J £ C^’‘^(M'^;M"^), i = 1,... ,m, j = and 7 E (1/3,1/2). Then 

we have 

\\X ’ — X ||-y,oo,[o,T] —^ 0 as n ^ 00 . 


This result directly implies: 

Corollary 1. Let a,bJ\c^^'^ E C'^’'^(M'^;M'^), i = j = 7 E (1/3,1/2) and 

G CI'\[0,T]-M.^). Then we have 

\\X ’ — X ||^,oo,[o,r] —^ 0 as n ^ 00 . 

Recalling that g^ = (id, VF, R^’”) Proposition [3] follows from Theorem [3] and the following two 
Lemmata. (Note that the following estimates are not covered by Section 15.5 in |3], since B^’^ is 
not a mollifier approximation.) 
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Lemma 3. For all 0 < 'j < < H there exists an almost surely finite random variable 

such that 

Proof. Clearly, it is sufficient to consider the one dimensional case. For t < 0 define = 0. Fix 
f, s € [0, T] and 7 ' G ( 7 , i7). First, consider the case |t — s| > Here, we have 



ft 


PS 

- Hf - B^’^ + Hf’^l < n 

/ (Hf - Hf )dn 

J t—lfn 

+ n 

/ (Hf - Hf )dn 

J s—l/n 


< 2\\B^ 




j <2 B^ 


l7',[o,r] 


Next, when If — s| < - one has 

’ I I — n 


\B^ -B^ - Hf’^1 < I + \B^’^ - 

< \\B^ '' " I ^ 


|t - sP' + n I / 


h',[o,T] 


< 2\\B^ 


<2 B^ vfnril^-^l 


lyjo,r] '^1 -"' 11 -^ ll7M0,T] I" '^1 „ 7 '- 7 - 


Lemma 4. Let 1/2 < 7 < 7 ' < Ff. Then, there exists an almost surely finite random variable 
.^t, 7 , 7 ' such that 


IG'^-GI 


C/([ 0 ,T] 2 ) 


< iFr,7,7' 


Proof. First, we prove the convergence the elements of the Levy areas which correspond to the 
smoothed fBm to the ones of fBm. Fix f, s G [0, T], i,j G {1, Clearly, we have 

J S J S 

< - Hf’«) - 

J S 

Set — Bu’'^'^''\ Applying the Young inequality with 1/2 < fj, < H yields 

[ (Z" - <Cf, \t-sf^, 

J S ^ ^ 

[ <Cf, \t-s\^^, 

J S ^ ^ 


i.e. we obtain 


t) 

\t — 


\ a a a a a a J 
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Now, we proceed with the parts that correspond to the iterated integrals which involve the 
Wiener process and the smoothed fBm. Fix i € {1,... ,m} and j € {1 ,... ,i}, s,t & Again, 

applying the Young inequality gives 


A(2)(s,t) = 

<Ga,m 


J s 




|t — s 


I 


with 0 < A < 1/2 ,0 < fj, < H such that A + /r > 1. It is only left to deal with 


A(3)(s,t) = 


J S J S 


But using the integration by parts formula for Young integrals we obtain 


f - W^) - j (W® - lFW)dZ, 

J s J s 

Using the previous step yields 


n,(j) 

u 


A(3)(s,t) <Ca,, 


Z'^Ao) 




_ oIA+a* 


|t — s 


The elements of the Levy area involving t and the “smoothed” fBms are easily treated. Here we 
have 


A(")(s,t) := 
A(5)(s,t) := 


f ^B^Ai) _ 

J S 

[ sd(Hf'« - 


< 

ZnAi) 

< 

ZnAi) 




|t — s 


|t — s 


1 + // 




Setting now 7 = (A + /r)/ 2 , the assertion follows from the previous lemma. 


□ 


5.2 Constructing numerical methods for mixed equations 

The almost sure convergence in the 7 -Holder norm of 

Yf" =xo+ [' (a(Yf+ c(Yf’")Hf’")ds + f 6 (Yf )d'W„ t € [0, T], (14) 

Jo Jo 

with 

= t€[ 0 ,r], 

to can be exploited to construct and to analyse numerical methods for mixed equations, proceed¬ 
ing similar to [mn]. In the latter references approximations of rough SDEs have been obtained by 
discretising their Wong-Zakai approximations. For example, applying an Euler discretisation with 
stepsize A = 1/n to (I14p yields the approximation 

Xk+i = Xk + a{xk)A + b{xk)(W(^k+i)A - + c{xk)iB^^ - H(^_i)a), k = 0,l,... 
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where = 0. Equation (|14l) is an Ito SDE with random coefficients, the only technical difficulty 
being the unboundedness of Using a localization procedure as e.g. in m and standard 

estimates involving the Ito isometry and Gronwall’s lemma one can show that 

I I P d.s. ^ 

sup l^kA ~ ^k\ —^ 0 , n^oo. 

k=0,...,\T/n] 

Corollary [1] then implies the convergence of this skewed Euler scheme, i.e. 

I hd I ^ d.s. 

sup \-^kA~^k\ —^ U, n —>■ oo. 

fc=0,...,rr/n] 

5.3 The natural Euler scheme for rough SDEs 

Using the correction formula, one can establish the convergence of the “natural” Euler scheme 

^ m 

Xk+1 = Xk+ (a{xk) + 2 X] ^ (15) 

i=l 

+ b{xk)(W(^k+i)A - WkA) + c(a:fc)(^^+i)A “ B^a), k = 0,l,... 

for the rough SDE 

= xo+ f aiX^)ds + f biX^)d^W, + f c{X^)d^B^, t G [ 0 , T], (16) 

JO Jo Jo 

at least for m = (. = 1. The notion “natural” is based on the following observations: for 5 = 0 
equation (fT 6 l) is an SDE driven by fractional Brownian motion with Hurst parameter H > 1/2, for 
which (fT5)l with 5 = 0 is a convergent approximation, see e.g. [miu, while for c = 0 equation (I16jl 
is a Stratonovich SDE, for which (|15l) with c = 0 is again a convergent scheme, see e.g. |11) . 

Using the results of |19| and Theorem [5] we have: 

Proposition 4. Let a,b,c € C'^’'^(M;M). Moreover let b'b G (^^’^(M;]^) and inf 3 ;gRc(x) > 0. Then 
there exists C > 0 such that 

sup {E\X^^ - Xkl^y^^ < C ■ 

k=0,...,\T/A] 


References 

[1] Kubilius, K. 2002. The existence and uniqueness of the solution of an integral equation driven 
by a p-semimartingale of special type. Stochastic Process and Applications 98(2): 289-315. 

[2] Mishura, Y., and Shevchenko, G. 2011. Stochastic differential equation involving Wiener process 
and fractional Brownian motion with Hurst index H > 1/2. Communications in Statistical 
Theory and Methods 40 (19-20): 3492-3508. 

[3] Lyons, T., and Qian, Z. 2002. System control and rough paths, Oxford University Press. 

[4] Friz, P., and Victoir, N. 2010. Multidimensional stochastic processes as rough paths. Theory 
and applications, Cambridge University Press. 


13 


[5] Mishura, Y., and Shevchenko, G. 2012. Mixed stochastic differential equations with long-range 
dependence: Existence, uniqueness and convergence of solutions. Computers and Mathematics 
with Applications 64(10): 3217-3227. 

[6] Shevchenko, G. 2014. Mixed stochastic delay differential equations. Probability Theory and 
Mathematical Statistics 2 (89): 181-195. 

[7] Goutin, L., and Qian, Z.2002 Stochastic rough path analysis and fractional Brownian motion. 
Probability Theory and Related Fields 122: 108-140. 

[8] Wong, E., and Zakai, M. 1965. On the convergence of ordinary integrals to stochastic integrals. 
Annals of Mathematical Statistics 36: 1560-1564. 

[9] Wong, E., and Zakai, M. 1965. On the relation between ordinary and stochastic differential 
equations. International Journal of Engineering Seience. 3: 213-229. 

[10] Karatzas, L, and Shreve, S. 1991. Brownian motion and stochastic calculus, 2nd edition. 
Springer, Berlin. 

[11] Kloeden, P., and Platen, E. 1999. Numerical Solution of Stochastic Differential equations, 3rd 
edition. Springer, Berlin. 

[12] Young, L. 1936. An inequality of the Holder type, connected with Stieltjes integration, Acta 
Mathematica by Uppsala 67: 251-282. 

[13] Gubinell, M. 2004. Controlling Rough Paths. Journal of Functional Analysis 216 (1): 86-140. 

[14] Friz, P., and Hairer, M. 2014. A course on rough paths. With an introduction to regularity 
structures, Springer, Berlin. 

[15] Deya, A., Neuenkirch, A., and Tindel, S. 2012. A Milstein-type scheme without Levy area 
terms for SDEs driven by fractional Brownian motion. Annales de I’Institut Henri Poincare. 
Probability and Statisties 48(2): 518-550. 

[16] Friz, P., and Riedel, S. 2012. Convergence rates for the full Gaussian rough paths. Annales de 
I’Institut Henri Poineare 50 (1): 154-194. 

[17] Jentzen, A., Kloeden, P., and Neuenkirch, A.2009. Pathwise approximation of stochastic dif¬ 
ferential equations on domains: Higher order convergence rates without global Lipschitz coef¬ 
ficients. Numerisehe Mathematik 112 (1): 41-64. 

[18] Davie, A. Differential equations driven by rough paths: An approach via discrete approxima¬ 
tion. 2008. Applied Mathematics Research eXpress , Article ID abm009, 40 p. 

[19] Mishura, Y., and Shevchenko, G. Rate of convergence of Euler approximations of solution 
to mixed stochastic differential equation involving Brownian motion and fractional Brownian 
motion.2011. Random Operators and Stochastic Equations 19 (4): 387-406. 


14 



